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The current bounds on the PPN parameters 7 and /3 are of the order of 10^^ — 10"^. Var- 
ious missions aimed at improving such hmits by several orders of magnitude have more 
or less recently been proposed like LATOR, ASTROD, BepiColombo and GAIA. They 
involve the use of various spacecraft, to be launched along interplanetary trajectories, 
for measuring the post-Newtonian effects induced by the solar gravity on the propa- 
, gation of electromagnetic waves. In this paper we investigate the requirements needed 

^y^' to measure the combination 1/ = (2 -|- 27 — /3)/3 entering the post-Newtonian Einstein 

, I ^ ' pericenter precession tj of a test particle to an accuracy of the order of 10~^ with a 

j^j^^ pair of drag-free spacecraft in the Earth's gravitational field. It turns out that the latest 

gravity models from the dedicated CHAMP and GRACE missions would allow to reduce 
the systematic error of gravitational origin just to this demanding level of accuracy. In 
■ regard to the non-gravitational errors, the spectral noise density of the drag-free sensors 

I^N , required amounts to 10^* — 10~^ cm s~^ Hz~2 over very low frequencies. Although 

^ ■ not yet obtainable with the present technologies, such level of compensation is much 

less demanding than those required for, e.g., LISA. As a by-product, an independent 
measurement of the post-Newtonian gravitomagnetic Lense-Thirring effect with a 1% 
accuracy would be possible as well. The forthcoming Earth gravity models from CHAMP 
and GRACE will further reduce the systematic bias of gravitational origin in both of 
such tests. 

Keywords: PPN parameters 7 and /3; drag-free spacecraft; CHAMP and GRACE Earth 
gravity models 



1. Introduction 

1.1. The PPN parameters 7 and (3 

In the Parameterized Post-Newtonian (PPN) formalisnff, developed by Nordtvedt 
and Will to make easier the comparison of metric theories of gravity with each 
other and with experiment, the Eddington-Robertson-Schiff parameters (3 and 7 
describe how much non-linearity is present in the superposition law for gravity and 
how much spatial curvature is produced by a unit mass, respectively; in general 
relativity [3 — ^ = 1. For a recent overview of the measurements of such PPN 
parameters see Refs. O [S] 



1 
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The most accurate and clean determinations of 7 come from the effects involv- 
ing the propagation of electromagnetic signals: in their PPN expressions only 7 
is present, at first order. For the Shapiro time delaj^, recently tested with the 
Doppler tracking of the Cassini spacecraft, the most accurate result i^ 7 — 1 = 
(2.1 ± 2.3) X 10~^, although the accuracy of such measurement has recently been 
somewhat questioned in Ref. [6] according to its authors, the motion of the Sun 
around the Solar System barycenter, if not properly accounted for, would induce a 
systematic bias of (S7 = 1.2 x 10"''. 

With f3 the situation is less favorable because it enters the post-Newtonian equa- 
tions of motion of test particles in conjunction with 7 itself in a rather complicated 
way. For example, in the post-Newtonian expression of the well known Einstein 
precessiorP of the pericentre w of a test particle the combination v = (2 -I- 27 — /3)/3 
is present. E.V. Pitjeva in Ref. [8] recently processed more than 317 000 Solar Sys- 
tem observations (1913-2003) for the construction of the Ephemerides of Planets and 
the Moon EPM2004 determining, among other things, j3 and 7 with an accuracy 
of 10^*. The same limit for /3 can also be obtained by com bining the Cassini result 
for 7 with the Lunar Laser Ranging tesP of the NordtvedlflS'liI effeclB, expressed 
in terms of 7y = 4/3 — 7 — 3: indeed, fronJ^^r; = 7 x 10"'', it follows SP = 1.8 x 10~^. 
Note that if = 1.2 x 10~^ is assumed^, the error in /3 beco mes 60 = 2.0 x 1 0"''. 

Various theoretical models involving scalar-tensor scenario a ***** ^ pre- 

dicts deviations from unity for f3 and 7 at a 10^^ — 10^^ level, so that it is very 

important to push the accuracy of such PPN tests towards this demanding accuracy. 

1 8 

The proposed LATOR mission ■ has, among its goals, the determination 
of 7 with an accuracy of 10^^ from the first order and a direct and inde- 
pendent measurement of /3 at a ~ 0.01% level via the second-order gravity- 
induced deflection of light. The expected accuracy in the proposed radio sci- 
ence experiments to be conducted with the future Mercury orbiter BepiColombo 
( http://sci.esa. int/science-e/www/area/index.cfm?fareaid=30) is of the order off^ 



10 ^ for /3 and 7. The accurate measurement of P and 7 is one of the scientific goals 
also of the proposed mission 

ASTROe20 21. the expected accurac y ia^ 1 x 10 
for hot]E . The astrometric mission GAIA should be able to reaclP^ ^7 ~ 10 
through the deflection of light. 



1.2. Aim of the paper 

A satisfactorily empirical corroboration of a pillar of physics like general relativ- 
ity requires that as many independent experiments as possible are conducted by 
different scientists in different laboratories. Now, all the performed or proposed 
high-precision tests of post-Newtonian gravity (in the weak-field and slow-motion 
approximation) are based on the effects induced by the gravitational field of the 



''It is a violation of the strong equivalence principle which is zero in general relativity. 
''See also Table 5 of Ref. |2] and Ref.l22l 
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Sun only. It is as if many independent experiments aimed to measure fundamental 
physical effects were conducted always in the same laboratory. Thus, it is worth- 
while to try to use different laboratories, i.e. other gravitational fields, to perform 
such tests, even if their outcomes should be less accurate than those conducted in 
the Sun's field. 

In this paper we will investigate if it is possible to realistically reach the 
10~'^ — 10^^ level of accuracy in measuring v with a couple of dedicated drag- 
free spacecraft in the terrestrial gravitational field. This analysis is motivated by 
the fact that, in regard to the systematic errors of gravitational origin, reaching a 
10~^ level would be right now possible thanks to the recent improvements in our 
knowledge of the classical part of the Earth's gravitational field due to the dedi- 
cated CHAMP (http://www.gfz-potsdam.de/pbl/op/champ/index_CHAM P.htmIl ) 
and GRACE ptp: // www.gfz-potsdam.de/pbl/op/grace/index_GRACE. htmll and 
http://www.csr.utexas.edu/grace/) missions. Much more demanding would be to 
reduce the non-gravitational errors down to the required level. As we will see, the 
drag-free technologies needed to implement such an ambitious goal are not yet 
available because they should work over observational time spans of several years; 
however, the y would be less demanding than those required for complex missions 
like LISaEI. 



2. The Einstein pericenter precession 

The Einstein pericenter precessiorP of a test particle in the gravitational field of a 
static central body is 

"inGM 
"'^^ - c2a(l-e2)' 

where G is the Newtonian constant of gravitation, M is the mass of the central body 
which acts as source of the gravitational field, c is the speed of light in vacuum, a 
and e are the semimajor axis and the eccentricity, respectively, of the test particle's 
orbit and n = ^jGM/a? is its Keplerian mean motion. 

As it is well known, Eq. |T]) allowed to explain the Mercury's anomalous per- 
ihelion advance of 42.98 arcseconds per century observed since the 19th century. 
Modern radar-ranging measurem ents of the planetary perihelia in the solar field 
yielded a 10"^ accuracji^SIMEll, 

Hiscock and Lindblom in Ref. [28] analyzed with some details the possibility of 
measuring the Einstein pericenter precessions of the natural satellites of Jupiter 
and Saturn which are much larger than those of the planets. lorio in Ref . [29] showed 
that, at present, this possibility is not yet viable, despite of the first data sets from 
the Cassini mission. 

The possibility of measuring the Einstein perigee rate in the Earth's gravity 
field with passive artificial satellites tracked via laser ranging was preliminarily 
investigated by Rubincam in Ref. [30] and Cugusi and Proverbio in Ref. [3TJ Ciufolini 
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and Ma tzner in Ref . [32l a nalyzed the data of the LAGEOS sateUite using the GEM- 
iPSEl and GEM T^M Earth gravity models claiming a total accuracy of 20%, 
but such estimate is unrealistic both because of the too optimistic evaluation of 
the impact of the classical part of the Earth gravity field and of the treatment of 
the non-gravitational forcea^Sj acting on the perigee of the LAGEOS-like satellites. 
Moreover, the very small eccentricity of the LAGEOS orbit (e = 0.0045) makes its 
perigee badly defined and highly affected by the non-conservative perturbations. A 
more accurate and reliable test was proposed in Ref. [37] by suitably combining the 
nodes of LAGEOS and LAGEOS II and the perigee of LAGEOS II, and using the 
EGM96 Earth's gravity modePSI. The estimated total accuracy is of the order of 
10"'^. An analysis of the use of the perigee of LAGEOS II can be found in Ref. [39l 
the total error is estimated to be 2%. Also the proposed LAGEOS-like satellite 
LARE^ could be used for measuring the Einstein perigee advancJ^. Hiscock 
and Lindblom in Ref. '28' envisaged the possibility of using geocentric drag-free 
spacecraft to measure the Einstein perigee precession. A step further concerning the 
use of dedicated terrestrial drag-free spacecraft for accurate tests of post-Newtonian 
gravity can be found in Ref. [42l In it the authors suggest that it would be possible 
to determine 7 and /3 independently of each other with an accuracy of 10"** by 
measuring v and rj. 

3. Outline of the various sources of systematic errors 

When a satellite-based test of a prediction of general relativity is proposed or per- 
formed, it is of the utmost importance to assess the total error budget in a very 
reliable and accurate way by accounting for all possible sources of systematic errors. 
Here we intend the aliasing action of other competing classical features of motion 
which may corrupt the genuine determination of the investigated relativistic effect. 
In regard to this topic, the choice of a suitable observable is very important in the 
sense that it must reach a good compromise between the errors of gravitational and 
non-gravitational origin whose simultaneous reduction is often in conflict. 

3.1. The gravitational perturbations 

The gravitational errors come from the mismodelling of the static and time- varying 
components of the multipolar expansion of the Newtonian gravitational potential 
of the Earth in spherical harmonicJ^. Indeed, the static part of the even (£ = 
2,4,6, ...) zonal (m = 0) harmonics J2, J4, Jg, ... induces secular precessionPH on a 
satellite's node and perigee iv 

^(gcopot) = ^l>2 ^-f-Jt^ 

(2) 

■^(goopot) = Z]^>2 

where ^Ij and Coj are the partial derivatives of the classical precessions with re- 
spect to the even zonal of degree I (see, e.g., Ref. l44| . while the odd zonal harmonics 
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J3, J5, ... affect the perigee with lonK-period harmonic perturbations having frequen- 
cies multiple of that of the perigee^. The secular variations of the even zon al har- 
nionics Ji, Jst ■■■ induce quadratic signatures on the node and the perigej ^^ ^ '^^ i 
There are also the long -period harmonic perturbations of the node and perigee 
due to the whose periods are determined by the lunisolar and satel- 

lite frequencies. The most insidious aliasing signals are the even zonal secular and 
quadratic advances and those harmonic signals that resemble superimposed linear 

trends over observational time spans of some years, like the I = 2,m = constituent 
4.7 

of the 18.6-year tide - : their mismodeled part can corrupt to a sign ificative extent 
the recovery of the genuine relativistic linear trends of interesiP ^^ ^ ^^ i 

A way to overcome the problem of the biasing classical even zonal secular trends 
consists in linearly combining the residuals (5^* of the Keplerian orbital elements, 
denoted generically as 5*, affected by the classical and relativistic linear rates 

s¥^) = ^ ¥fdJ, + vi/«e„i, I = N, (3) 

where N is the number of perigees and/or nodes employed, in order to solve for the 
parameter ^j-ci, which accounts for the relativistic effect of interest (^roi = in our 
case) , and cancel out the impact of as many e ven zonal harmonics as possible among 
those recognized to be the most biasing oneJ^^EHSl]^ rpj^g coefficients ^^^-^ are the 
partial derivatives of the classical precessions of with respect to the even zonal 
of degree £. The so-obtained linear combinations are sensitive to general relativity, 
irrespectively of the level of mismodelling dJf, of the static and time-varying parts 
of the selected even zonal harmonics. It turns out that the accuracy reached in 
the latest Earth gravity models from CHAMP and GRACE would allow right now 
to reduce the systematic error due to the geopotential to 10~^ of the Einstein 
precession if a suitable orbital geometry is chosen for a dedicated spacecraft-based 
mission. We discuss this topic in detail in Section [H 



3.2. The non- gravitational perturbations 

On the other hand, such really important benefits cannot be fully capitalized if the 
non-g ravitational forced acting on the spacecraft like the direct solar radiation pres- 
sure, Earth's albedo. Earth's IR radiation, atmospheric drag, thermal effects, etc^^Hl 
are not suitably accounted for and reduced to the same level of the gravitational 
bias. For example, in the case of the Lense-Thirr ing^ tes ts with the existing LA- 
GEOS and the proposed LARES/OPTIS satelliteP^iS4i55|^ ^^le obtainable accuracy 
is set up to ^ 1% by the non-conservative forces acting on the passive LAGEOS and 
LAGEOS II satellites which are not endowed with a ny active mechanism of compen- 
sation of such kind of perturbations [56| 57 | 58 | 59 [601 ^j^g systematic error due to the 



'^Let us recall that their magnitude is proportional to the area-to-mass ratio (r/m of the satellite and 
depends on the interaction of the electromagnetic radiation with the satellite's physical structure, 
on its rotational state and its orbital geometry often in a very complex and intricate way. 
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geopotential would be much smaller according to the latest CHAMP/GRACE-based 
modelJ^. Thus, it is apparent that only drag-free spacecraft could push the bias 
due to the non-conservative forces to the same level of that due to the gravitational 
perturbations. This topic will be the subject of Section [6l 



4. How to reduce the gravitational perturbations 

The idea of using a pair of drag-free Earth's artificial satellites was considered also 
in Refs. [SH [531 IMl [Ml in that case the investigated feature of motion was the Lense- 
Thirring force which affects the node and the perigee of a satellite (see Section [5]) . 
In Refs. [62[ [63l [64l it was shown that the difference of the perigees of two satellites in 
supplementary orbital configuratior|^ is able to add the Lense-Thirring precessions 
and to cancel out the classical even zonal precessions because of their functional 
dependence on the inclinatiorPSI. The use of active mechanisms of compensation 
of the non-gravitational perturbations was suggested because of the strong impact 
that they have on the perigees of LAGEOS-like satellites. 

In the case of the Einstein precessions the situation is different since they are 
independent of the inclination. This implies that the difference of the perigees in 
the supplementary orbital configuration is not a good observable because it would 
exactly cancel the relativistic rates as well, independently of the values of 7 and (3. 
The linear combination approach involving the nodes and the perigees of both the 
spacecraft must, thus, be followed. Having the residuals of four measurable Keple- 
rian orbital elements at our disposal it is possible to design a suitable observable 



•,(1) 



■ (2) 



(4) 



which is affected by the Einstein precessions and is independent of the static and 
time-dependent parts of the first two even zonal harmonics of geopotential J2, ^4, 
and of the Lense-Thirring precessions. It turns out that, also in this case, the sup- 
plementary orbital configuration is not suitable because the coefficients weighing 
the perigees would be equal and opposite cancelling the relativistic signature. If 
different orbital parameters are chosen for the two spacecraft the situation be- 
comes favorable. Indeed, with the orbital parameters of Table [H we have for our 



Table 1. Orbital parameters of the proposed drag- free 
spacecraft (1) and (2): a is in km and i in deg. Also 
their relativistic precessions, in mas yr~^, are listed. 



Satellite 


a 


i 


e 








(1) 


13000 


50 


0.2 


2955.7 


-53.0 


27.5 


(2) 


12000 


103 


0.2 


3610.5 


23.6 


35.0 



"^It means ai = 02, ei = £2, U = *2 — 180 deg. See Refs. [65l 1661 
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J2 — J4— Lense-Thirring-free combination 
' Cl = 2.2817, 

C2 = -0.6861, 

(5) 

C3 = 0.5165, 

^ u'^l + ciuj^qI = 11194.1 mas yr"\ 

and the relative systematic error due to the uncanceled even zonal harmonics 
amounts to 

- ^ 3 X 10-^ (6) 

^ gcopot 

if the currently available CHAMP/GRACE-based EIGEN-CG03C Earth gravity 
model^is adopted for the calculation. Note that Eq. ([6]) represents a realistic upper 
bound of the gravitational bias because it has been obtained by linearly summing 
the absolute values of the individual errors up to degree i = 20. Note that this 
level of truncation is quite adequate because of the high altitudes of the satellites: 
indeed, it turns out that the numerical value of ^'^/'^Igcopot does not change if one 
adds more even zonals. The current models from CHAMP and GRACE are still 
preliminary and do not yet exploit the full data sets gathered by such missions 
which are still ongoing. Thus, it is likely that when new, more robust and reliable 
Earth gravity solutions will be available the estimate of Eq. ([6|) will further improve, 
all the satellites' orbital parameters being equal. 

In regard to the time-varying gravitational errors, the periods of the nodes and 
the perigees of the satellites are listed in Table [2j It can be noted that they amount 

Table 2. Periods, in yr, 
of the node and the 
perigee of the proposed 
drag-free spacecraft (1) 
and (2) 



SatelHte 


P(0) 




(1) 


-1.71 


2.06 


(2) 


3.7 


-2.22 



to just a few years. This means that the sinusoidal perturbations induced on the 
nodes and the perigees by the ^ = 2,m = 1 constituent of the Ki solar tidJ^ 
and the odd zonal harmonic J3, which have just the same periods of the satellite's 
node and perigee, respectively, and are not canceled out by the combination of Eq. 
(j4|), do not represent serious bias over reasonably long observational time spans. 
A potentially less favorable situation occurs for the £ = 3,m — l,p ~ l,q = — 1 
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oceanic constituent of the Ki tidJ^ which is not canceled by Eq. ^ and affects 
the perigees with perturbations whose periods are -10.0 yr and -5.6 yr, respectively. 
For the complete details of the action of the £ — 3, m — 1 oceanic constituents of 
Ki on the perigees of (1) and (2), see Table [3l According to the global ocean tide 

Table 3. Periods P, in yr, and mismodeled amplitudes SA, in mas, of the tidal per- 
turbations induced on the perigees of the proposed drag-free spacecraft (1) and (2) by 
the £ = 3, m = I oceanic constituents of the Ki solar tide according to the results of 
Table 6.4.8.5-1 of Ref. 36. 



Satellite 


P{p=lq = -1) 


SA{p = 1 g = -1) 


P(p = 2 <7 = 1) 


5A{p = 2 g = l) 


(1) 


-10.0 


-10.1 


-0.9 


2.1 


(2) 


-5.6 


7.0 


1.4 


3.3 



model of the EGM96 solutiorPSI^ it turns out that the most insidious tidal effect 
is the Ki £ — 3, m — 1, p = 1, q — —1 perturbation on the perigee of satellite 
(1) with a mismodeled amplitude of -10.1 mas. Thus, the maximum relative error, 
over 10 yr, on the combination of Eq. ([4]) amounts to 9 x 10~^. However, it must be 
noted that such an estimate is largely pessimistic both because it assumes a peak 
value after 10 yr and because it is based on the relatively old EGM96 model: more 
recent and accurate models like GOT99, GSR 4.0, FES2004 would greatly reduce 
the impact of such a source of bias. 

If in constructing the combination of Eq. ^ one chose to cancel Jg instead of the 
Lense-Thirring precessions, and a nearly-polar orbital configuration was adopted 
for both the satellites, a bound iJi^/i^lgcopot — 10"'' could be reached, according 
to EIGEN-CG03G. But, in this case, it would not be possible to separate^ the 
Einstein precessions from the Lense-Thirring ones which are weighted by a different 
combination of 7 and B, namely fi = (1 + 7)/2. Moreover, inclinations close to 90 
deg would be harmful'"' because of the extremely long node periods and, thus, of 
the Ki tidal perturbations. 

5. The impact of the inclination errors 

Another source of systematic error is related to the uncertainty in the knowledge of 
the spacecraft inclinations and their impact on the bias due to the uncanceled even 
zonal harmonics: indeed, recall that the coefficients and of the classical even 
zonal precessionJ331 do depend on the inclinations as well. 

To this aim, we considered ^'^/'^Igoopot ^ function of the two variables i^^'' and 
i^^^ and studied its behavior in small ranges Ai^^^ and Ai*-^-* around the nominal 
values of Table [1] which simulate the errors in the inclinations. It turns out that for 

"^Also the post-Newtonian J^lf? effects due to the Earth 's oblatenesJSHIMl affect all the combi- 
nation proposed here. They must be modeled, otherwise they would induce a systematic bias of 
10~^ on the measurement of the Einstein precession. 
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a very conservative and pessimistic assumption Ai = 2 deg for both z^^^ and i'-^-' 
the relative error in v varies by just 9 x 10~^. 

6. The accuracy required for the non-gravitational perturbations 

hi this Section we will investigate the level to which the non-gravitational accelera- 
tions should at most affect the spacecraft if a relative systematic error of 10^® was 
to be obtained. 

The Gauss equations for the variations of the node and the perigee arJ 



<Kl _ 1 

duj dfl vT" 

— — — coai— — I 

at at m 



A. 



„ (^) sin(c. -t- /), (7) 



-Ar COS / + ( 1 + - ) sin / 



r 



(8) 



in which / is the true anomaly counted from the pericentre, p = a(l — e^) is the 
semilactus rectum of the Keplerian ellipse, i is the inclination of the orbit to the 
Earth's equator, A^, At, An are the radial, transverse (in-plane components) and 
the normal (out-of-plane component) projections of the perturbing acceleration A, 
respectively, on the co-moving frame {r,t,n}. 

We are interested in the effects averaged over one orbital period P: they can be 
obtained by evaluating the right-hand-sides of Eq. ((T]) and Eq. ([8]) on the unper- 
turbed Keplerian ellipse 

V 



and averaging them with 



1 + e cos / ' 
dt _ (1 - e2)3/2 



(9) 



■df. (10) 



P 27r(l -|-ecos/)2 

Let us start with the node and consider the action of an out-of-plane acceleration 
consisting of a part TV'"' constant over one orbital revolution and a part varying 
with the orbital frequency 

A„-Af^"^+Afssin/ + iVccos/, (11) 

where Ng and N(j are to be considered as constant over one orbital revolution. Now 
Eq. (jlip must be inserted into Eq. ([7|) and treated as previously outlined. By using 

/ r2TT _sin(uj+f),r _ 37re • 

Jo (l+ecos/)a"j ^ (l-e2)5/2 Wnw, 

< 

. Jo (l-Hecos/)3 - (l-e2)5/2 [i-L " 6 jiVs COS w + (L + 2e )l\c sm 

one finally gets 

(f ) = jkl^cmh^i [(1 + 2e^)^c - 3eiV(0)] sin^ + (1 - e')Ns cos^} . (13) 
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The same procedure holds for the perigee. By inserting 
Ar = + Rs sin f + Rc cos /, 

At = T(") + T5 sin Z + Tc cos/ 
into the second term of the right-hand-side of Eq. ^ and using 

/ r27r COS / u _ 2ire 
Jo (T+ecosTF (l-e2)3/2i 

.2. (2+e cos/) sin/ ^^ ^ q 
JO (l+ecos/)'=' ' 



(14) 



r2?r cos /(-Rg sin f+Rc cos /) 



(T+ecos/p 



df 



27tRc 



(l_e2)3/2 



(15) 



/c 



271- (2+e 



5 /) (Ts sin /+Tc cos/) sin / ^^ _ 
(l+ecos/)i' ".^ 



it is possible to obtain 



= — cos I 



dn 



dt I 



dio \ 



(l_e2)3/2 



(C/S) 



ith 



v(0) 



/ \ 
\ dt / 

/dw\(C/S') 
\ dt / 



a(l-e2) 
GM 



(16) 



(17) 



1 - 



Rc 



1 



2(l-e2)3/2 



(l_e2)3/2 

This result tells us that the perigee would be affected by secular or long-period 
perturbing rates, according to the frequencies embedded in N'-^\ Ns, Nt, R'-°\ Rc 
and Ts- 

From Eq. (fT7|) it is possible to obtain the level to which the non-gravitational 
perturbations should be reduced in order not to exceed a given systematic error in 
the measurement of the gravitoelectric perigee precession. For a ~ 12000—13000 km 
and e = 0.2, a 10~^ relative accuracy could be reached with a residual acceleration 
level o^ 10^^^ cm s~^. The drag-free technologies currently under development 



for, e.g., the LISA mission^ should allow a cancelation of the non-gravitational 
accelerations down to 3 x 10""'^^ cm s~^ Hz~ 2 . In our case, with an orbital frequency 
of 7 X 10~^ Hz, a drag-free level of ^ 10~^ cm s~^ Hz~2 would be sufhcient: such 
a cancelation should be reached and maintained over time spans some years long, 
i.e. in a frequency range spanning from about to 10^^ Hz. It is certainly not an 
easy task to reach. No long-term experience about the performance of such systems 
is available at present; current high pre cisio n accelerometers show-by extrapolating 
present data-a spectral noise density offl^ 10~^ cm s^^ Hz^^. A very accurate 
in-orbit calibration would be required during the entire duration of the mission. 



*For a comparison, the impact of the direct solar radiation pressure on the existing LAGEOS 



satelhtes amounts ti 



10" 
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7. The measurement errors 



For the perigee and the node the measurable quantities are r — eauj and r = a siniQ. 
Thus, for the combination of Eq. ^ the measurement error can be written af|f| 



1 



1 



1 



1 



a(i)e(i) 



■ ci 



a(2)e(2) 



C2- 



(1): 



ea- 



rn 



a(2) sini(2) 

By using the values of Table [T] and by assuming (5ri„cas = 0.1 cm over an observa- 
tional time span T = 10 yr, one gets for the relative error 



= 2 X 10" 



(19) 



About the assumptions made, it must be noted that we consider an experimen- 
tal millimeter accuracy in determining the spacecraft orbit, in a root-mean-square 
sense, over 10 years: it is certainly not an easy task to be accomplished, especially 
in view of the fact that no empirical along-track accelerations should be fitted: oth- 
erwise the signal of interest could be removed as well. Moreover, it must be said 
that attempts to improve measurement errors are not the same as the capability of 
reconstructing the orbit with mm accuracy over 10 years. 



8. The measurement of the Lense-Thirring effect 



As a by-product, our active spacecraft could a lso be used to make an independent 
measurement of the Lense-Thirring 

precessioniSSESESIIllia 



LT 



'^LT = -M 



„ 2&S 

^^c2a3(l-e2)3/2 I 

6GS cos i 



(20) 



c2a3(l-e2)3/2 ' 

where S is the spin angular momentum of the central rotating mass. Indeed, if 
one combines the nodes and the perigees so to cancel out J2, and the Einstein 
precessions, and solves for /i, the resulting combination is obtained 



LT 



with 



; • (2) 



2"lT 



^(2)1 



3"ltJ 



(21) 





-0.9, 


k2 = 


0.21, 




= -2.7, 



(22) 



,(1) 

■^LT 



; • (2) 



^LT 



S(2) 
LT 



T61 mas yr 



Please note that the labels (1) and (2) refer here to the two spacecraft. 
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In this case, the major hmiting factor is represented by the systematic error due 
to the uncanceled even zonal harmonics because the drag-free level of cancelation 
of the non-gravitational perturbations discussed in Section [6] would allow to reduce 
their impact well below the percent level. It turns out that (5/x//L(|geopot = 9 x 10~^, 
according to EIGEN-CG03C. Such a result is not so good as it would be with 
the supplementary orbital configuration, and it could be reached also with only 
one new passive satellite and the existing LAGEOS and LAGEOS II (see Ref. I6T|) . 
Nonetheless, (5/i//i|goopot should be improved by the new Earth gravity field solutions 
from CHAMP and GRACE pushing the overall accuracy of such test towards the 
~ 0.1% level. This achievement could not be reached with the low-cost alternative 
involving one new satellite and the LAGEOS ones because of the uncanceled action 
of the non-gravitational perturbations. 

9. Conclusions 

In this paper we investigated what is needed to measure the combination (2 + 
27 — /3)/3 of the post-Newtonian gravitoelectric Einstein perigee precession of a 
test particle at a « 10~^ level of accuracy by using a couple of drag-free spacecraft 
in the Earth's gravitational field. The observable proposed is a linear combination 
of the nodes 57 and the perigees lo of such satellites to be analyzed over an ob- 
servational time span of several years. It allows to cancel out the impact of the 
static and time- varying parts of the first two even zonal harmonics of the terrestrial 
geopotential, which represent a major source of systematic error, and of the post- 
Newtonian gravitomagnetic Lense-Thirring precessions. By using the latest Earth 
gravity model EIGEN-CG03C from CHAMP and GRACE, the bias induced by the 
remaining uncanceled even zonals amounts to 10~^, for a suitable choice of the or- 
bital configuration of the satellites with a ^ 13000 — 12000 km and i ^ ~ 103 
deg. When the forthcoming, more accurate and reliable CHAMP/GRACE-based 
Earth's gravity models will be available, it should be possible to push the gravita- 
tional bias down to 10~^. An active reduction of the impact of the non-conservative 
perturbations is required in order to fully benefit of the low systematic error of grav- 
itational origin. It turns out that the level of cancelation of the non-gravitational 
perturbations should be of the order of 10~^ cm s~^ Hz~2 . Although the currently 
available drag-free sensors do not yet allow to reach such performances, also because 
they should be maintained over several years, they are less demanding than those 
required for, e.g., LISA. As a by-product of the proposed mission, an independent 
measurement of the Lense-Thirring effect with a « 1% accuracy would be possible 
as well. In this case, given the drag-free level required for the previous test, the 
most limiting factor is represented by the uncanceled even zonals; the forthcoming 
CHAMP/GRACE models should reduce their impact on this test as well. Another 
stringent and demanding requirement is reaching the ability of reconstructing the 
satellites' orbits at wmm level in root-mean-square sense over a time span of 10 
years. 
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